The role of surface plasmon polaritons in the optical response of a hole pair 
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The optical emittance of a hole pair perforated in an opaque metal film is studied from first- 
principles using the coupled-mode method. The geometrical simplicity of this system helps to un- 
derstand the fundamental role played by surface plasmon polaritons (SPPs) in its optical response. 
A SPP interference model without fitting parameters is developed from the rigorous solution of 
Maxwell's equations. The calculations show that the interference pattern of the hole pair is deter- 
mined by two scattering mechanisms: (i) re-illumination of the holes by the in-plane SPP radiation 
and (ii) an effective impedance depending on the single-hole response. The conditions for construc- 
tive and destructive interference only depend on the phase difference provided by each of the two 
scattering mechanisms. 



I. INTRODUCTION 

The extraordinary transmission through nanohole ar- 
rays milled into metallic films [l[ is attributed to the res- 
onant excitation of surface electromagnetic (EM) modes 
by the incident light Q. In the optical regime, these sur- 
face EM modes are surface plasmon polaritons (SPPs), 
modified by the metal corrugation. The light-SPP cou- 
pling is made possible by the additional grating momen- 
tum provided by the scattering of the incident light by 
the hole array. Nevertheless, interference of excited SPPs 
is set up even for two interacting holes 0t3]- Increasing 
the number of holes, the transmission is enhanced due to 
better-defined peaks of the structure factor, appearing at 
the reciprocal lattice vectors Q. It must also be noted 
that light-SPP interaction is not the single mechanism 
behind the extraordinary optical transmission (EOT). 
The EOT physical scenario is completed by the excita- 
tion of localized and Fabry- Perot modes J9rll61 ]. which 
may also contribute the whole process (see [17| for a com- 
prehensive review). 

The aim of the present paper is to study the interfer- 
ence pattern of the simplest interacting system: a hole 
pair. Since the original proposition of the "nanogolf" 
effect by Sonninchscn et al. [3| , several groups have mea- 
sured the optical interaction of two holes, see for exam- 
ple [g, @, [l8( ■ These groups have used basic SPP reso- 
nant models in order to explain the characteristic optical 
transmittance of the hole dimer, which oscillates as func- 
tion of the hole-hole distance, with period equal to the 
SPP wavelength. These approaches have in common that 
the relevant scattering channels are assumed ad hoc: only 
SPP scattering channels are included in the final optical 
response. 

In this paper we make no such assumption and solve 
Maxwell's equations from first-principles using a coupled- 
mode method (CMM) [U Gj|. We shall consider the 
two possible radiative channels: freely propagating light 
radiated out-of-plane into the far field, and SPP power 
scattered along the metal plane. The out-of-plane power 



FIG. I. (Color). Schematic representation of the hole-pair 
geometry. 



P ra d, normalized to the power incident on the hole area, 
gives the far field transmittance T. This is the quan- 
tity commonly used to characterize EOT. However, to 
the best of our knowledge, the in-plane SPP power P spp 
has not yet been measured for a hole pair. We shall ana- 
lyze the relevant scattering mechanisms for each radiative 
channel. Moreover, we shall derive, without fitting pa- 
rameters, the conditions for constructive and destructive 
interference, hereafter conditions for interference (CI) , 
that explain experimental interference patterns [a, LL8| . 

The paper is organized as follows. In the next section 
we briefly review the CMM and give the expressions for 
P spp and Prad- The assumptions behind the CMM and 
some cumbersome mathematical formula are reported in 
the Appendix. For the sake of completeness, section IIIII 
summaries the emittance of a single hole. Section IIVI 
discusses the optical response of the hole pair. A sub- 
section is devoted to clarify the scattering mechanisms 
dominating the conditions for interference. At the end, 
we outline the main conclusions of the paper. 



II. THEORETICAL FRAMEWORK 

Figured] renders the hole-pair geometry studied in this 
paper. Two identical circular holes of radius r^, sepa- 
rated by a distance R, are milled into an infinite metal 
film of thickness h and dielectric function e m . In general, 



the metal film lays on a substrate with dielectric con- 
stant e, it is covered with a dielectric superstate ei, and 
the space inside the holes is characterized by a dielectric 
constant e 2 . For the sake of simplicity, e = ei = 62 = 1 
is used along this paper. We consider in what follows 
that the metal film is illuminated by a normal-incident 
p-polarizcd plane wave, oriented along the main axis of 
the hole pair, as shown in Fig. [T] We shall focus on 
the energy power radiated into the transmission region 
(z>0). ' 

Maxwell's equations are solved self-consistentl y u sing 
a convenient representation for the EM fields PJI, 13 ■ In 
both substrate and superstrate the fields are expanded 
into an infinite set of plane waves with both p- and s- 
polarizations. Inside the holes the most natural basis is 
a set of circular waveguide modes. Convergence is fast 
achieved with a small number of such modes [2CJ, [21] . In 
fact, we shall see that the fundamental waveguide mode 
is a good approximation for our problem. The assump- 
tions behind this coupled-mode method, as well as its 
relevant constitutive quantities, are briefly review in the 
Appendix under the single mode approximation. 

The flux power traversing the hole is distributed into 
two channels [131 : W out-of-plane radiation, freely prop- 
agating into the far-field, and (ii) SPP power, scattered 
along the metal plane. The calculation of these two quan- 
tities is straightforward within the CMM after we know 
the amplitude of the fundamental waveguide mode at the 
hole openings E[, where i = 1,2 labels each hole. For a 
normal-incident plane wave, both holes receive the same 
illumination / (|A1|) . therefore E[ = E' 2 = E' due to the 
symmetry of the system with respect to the central point 
of the hole pair; E' hence reads 
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(1) 



where the hole- hole propagator Ghh{R) (|A4[) represents 
the coupling of the two holes as a function of the hole- 
hole distance R. This interaction can be seen as a re- 
illumination of the hole i by the magnetic field GhhEj 
radiated from the other hole j. Notice that there is also 
a self- illumination term for each hole, G s h, which adds 
to the single-hole scattering mechanisms E (|A2[) and G v 
(|A3[) . Using Eq. ([T]), the out-of-plane power emitted by 
the hole pair simplifies to 



P rad (i?) = |£'| 2 <7 rad (i?), 



(2) 



sh 
Srad 



.CdW Provides 



where the propagator g la d(R) 

the far field radiated from the hole pair, gf^ represents 
the contribution of each single hole, and g"l^(R) (|A8[) is 
a term arising from the interference of the fields radiated 
by the two holes. 

On the other hand, we can obtain the power radiated 
into SPPs by computing the contribution from the plas- 
mon pole in the propagator [l9| ■ The power of the scat- 
tered SPPs is first computed, at a point r on the metal 
surface several SPP wavelengths away from the nearest 
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FIG. 2. Normalized-to-hole-area out-of-plane (P ra d) and in- 
plane SPP (P sp p) emittance as function of the hole radius rh 
(in nm) for a single hole milled in a silver film, free standing 
on air (h = 250 nm and A = 700 nm). Symbols and lines 
represent converged results and the single mode approxima- 
tion, respectively. The inset show the cutoff wavelength, A c 
(in nm). 



edge of the hole pair, by integrating the in-plane ra- 
dial component of the Poynting vector, defined with the 
SPP fields, on a cylindrical surface of radius r and semi- 
infinity extension in z > 0; the power in the plasmon 
wave is then calculated using the known decay length of 
the SPP. The integrated power reads 



P SPP (R) = \E'\ 2 g spp (R), 



(3) 



where the propagator g spp (R) = g B spp + ff^(-R) provides 
the total SPP field radially scattered along all possible 
angular directions in the metal plane, g s ^ p (|A11[) repre- 
sents the contribution of each single hole, and g" pp (R) 
(|A12|) is the interference term 

The conservation of the energy flux for a lossless metal 
(|A15[) . imposes a constrain to the real part of the full 
interaction propagator G(R) = G s h + Ghh(R), which ful- 
fills Re[G(R)] = g Ta d + g spP , see Appendix. We shall see 
that this relation is a good approximation for lossy met- 
als at optical frequencies, and we shall use it in section 
IIVI However, the hole-hole interaction can not be fully 
understood without a previous knowledge of the optical 
response of a single hole, which is briefly reviewed in the 
next section. 



III. SINGLE HOLE EMITTANCE 

The emittance spectrum of a single circular hole is de- 
scribed in this section for the sake of completeness, al- 
thou gh t his issue have been largely study, see for example 
H S Hi 0, EE EHH and references therein. The be- 
havior of both P ra( j and P spp is depicted in Fig. [5] as 



function of the hole radius, for a free-standing Ag film 
with h = 250 nm. The Ag dielectric function e m (A) is fit- 
ted to Palik's data [231. ^ is equal to e m = —19.9 + 1.15 i 
for the incident wavelength, A = 700 nm, which is kept 
constant along the paper. Both P ra( j and P spp are nor- 
malized to the power incident on the hole area. 

Fig. [5] renders P spp for both the fundamental mode 
approximation (dashed line) and converged results (open 
circles). Both curves practically overlap, so the funda- 
mental mode is enough to achieve converged results for 
this emittance channel. For the out-of-plane emittance 
the agreement between single mode (solid line) and full 
calculations (full circles) is slightly worse, but still the 
difference is less than 15% and tendencies are well cap- 
tured in the parameter range considered. 

As already stressed in Ref. [2CJ, P spp irh) presents a 
broad peak with maximum at rh = 190 nm, close to the 
cutoff radius, r c = 168 nm for A = 700 nm. The cutoff 
wavelength, A c , is represented in the inset of Fig. [2~] as a 
function of r^ . The resonance appears in the field at the 
opening \E'\ (not shown), while the decay for r^ > r c is 
due to that in the single hole SPP propagator gfl (|A11I) . 
For r^ > ?' c most of the the energy is radiated out of the 
plane. In this case, both gf^ d (not shown) and P ra d reach 
a fast saturation with the hole radius. 



IV. OPTICAL RESPONSE OF A HOLE PAIR 

We define the normalized hole pair emittance as the 
power radiated into each channel, out-of-planc ([2]) or in- 
plane-SPP (J3j> , divided by twice the corresponding emit- 
ted power of a single hole located at R = 0, i.e. 
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where E' N is the ratio of the electric field at the hole 
openings 
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and we have used that the illumination of an isolated 
hole is equal to the illumination of each hole in the pair 
for a normal incident plane wave; while the ratios for the 
out-of-plane and SPP propagators are given by 
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The normalized omittances i~j ra d and rj spp are depicted 
in Fig. |3] as a function of the hole-hole distance R for 
increasing hole radius; r^ = 100 nm (blue dashed line), 
150 nm (red solid line), and 250 nm (black short-dashed 
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FIG. 3. (Color). Normalized out-of-plane emittance r\ ra d (a) 
and normalized in-plane SPP emittance r\ avv (b) as function of 
R/Xspp for increasing rh', rh = 100 nm (blue dashed line), 150 
nm (red solid line), and 250 nm (black short-dashed line). The 
holes are milled in a free standing Ag film of thickness=250 
nm. The illumination wavelength is A=700 nm. The CI re- 
ported in IIV Al are included in both (a) and (b) for r^ — 150 
nm: maxima (at R/\ spp = m — 1/4 for r\ ra d and R/\ spp ~ m 
for n S pp) and minima (at R/\ spP = m + 1/4 for rj ra d and 



R/Xs 



1/2 for n spp ) are represented with arrows, 



while vertical dashed lines are used for the condition r\ = 1 
(at R/Xsjjp — (m + l)/2 for rj rad and R/\ 3VV ~ (2m + l)/4 
for rj spP ); m — 1,2,3,... 



line). R is normalized to the SPP wavelength X spp = 
2ir/Re[k spp ], where k spp (|A10[) is the SPP propagation 
constant in silver; \ spp — 682.3 nm for the chosen A = 
700 nm. 

In accordance with experimental works [6j, l7J , the com- 
puted powers ri ra d and rj spp oscillate with period X spp . 
However, rj ra d behaves different than ij spp as a function 
of the hole radius. The amplitude of rj ra d strongly os- 
cillates with r/j. Indeed, increasing r^ from 100 nm to 
150 nm at fixed R we can transform a maximum of r\ ra d 
into a minimum. To the best our knowledge this depen- 



dence of r\ ra d on r^ has not been previously reported. 
Moreover, in the thin-film limit it has been found that 
the CI only depend on the edge-edge distance, and not 
on r/i . Further experimental work is needed to study the 
dependence on r^ for opaque metal films and hole sizes 
larger that the metal skin depth (the region of the pa- 
rameter space targeted in this paper). Nevertheless, it 
is worth stressing that the available experimental data 
IS 01 report the same CI for very different systems [2J] . 
In both cases r^ is very small (~ A/20), but while Ref. 
@ considers a thin gold layer (h = 20 nm) on a glass 
substrate, Ref. Q uses an optically thick silver film, 
immersed in a medium with refractive index n = 1.45. 
Both experimental CI are the same as for a third dif- 
ferent system, the particular case r% = 150 nm in Fig. 
[3] (a), i.e. maxima arc at R/X spp = m — 1/4, minima 
at R/X spp = m + 1/4 (both represented with arrows), 
and r\ ra d = 1 at R/X spp — (m + l)/2 (represented with 
vertical dashed lines), where m = 1, 2, 3, . . . 

In contrast, the amplitude of r] spp shows a stronger de- 
pendence on R, but does not present such large variations 
with size of the holes. Maxima of rj spp occur close to the 
conditions for constructive interference of SPPs at the 
flat metal surface (R = TnX spp ), while minima appear 
close to conditions for destructive interference of SPPs 
between the holes (R = (2m — l)X spp /2). As the energy 
traversing the holes is distributed into the out-of-plane 
and in-plane channels (|A15[) . we find in Fig. [3] that r] ra d 
and r] spp behave as complementary scattering channels, 
with a relative contribution that changes as a function 
of both R and tv Taking into account the interference 
pattern on the SPP channel we solve the apparent para- 
dox put forward in Ref. [7(: although r\ ra d is described 
by an SPP interference model, there is neither a trans- 
mission enhancement nor suppression at the conditions 
for constructive interference of SPPs. 

The out-of-plane radiation, r\ ra d (|4]), is mainly deter- 
mined by E' N , i.e. by the change in the field at the hole 
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FIG. 4. (Color), (a) Normalized out-of-plane emittance rjrad 
(red solid line), its constituent terms |i?iv| 2 (black short- 
dashed line) and <? rad (dark- yellow dashed line), and the ap- 
proximate expression for rjrad (orange dash-dotted line) of 
Eq. p0|) . (b) Normalized SPP emittance r) spp (blue solid 



due to the presence of the other hole. This is illustrated line )> ,j ts constituent terms \E' N \ (black short-dashed line) 

and g^pp (magenta dashed line), and the approximate expres- 
sion for rj S pp (violet dash-dotted line) of Eq. I)13p . All these 



in Fig. Eta), where r) ra< i is compared with both l-E^I 2 and 
■9rad- We observe that the interference between the radia- 
tive field of the two holes, given by gf ad , practically does 
not change the total transmission for hole-hole distances 
larger than 2X spp . Conversely, the normalized in-plane 
propagator gf plays an important role setting up the 
CI for r] spp ((5]), see Fig. HJb). Although the contribution 
of |-E^v| 2 can not be neglected, the interference pattern 
of g^ pp resembles the behavior of rj spp . 

The CI developed in the next section strongly depend 
on the properties of the in-plane propagator Ghh(R), 
which is behind the interference pattern of both radia- 
tive channels. We use the following decomposition of the 
in-plane propagator 



quantities are represented as function R/Xspp- The hole ra- 



dius is r/j 
in Fig. H 



150 nm, the rest of parameters are the same as 



Ghh(R) 
where G^ d (R) 



Ghh 
rad 



(R) 



G^JR) + G h c !:(R), (9) 



represents the contribution of ra- 
diative modes, G^ V (R) (|A9[) designate the contribution 



spp \ 



of the plasmon pole to evanescent modes, and G^(R) 



(|A13[) denotes the contribution of the remaining evanes- 
cent modes. This decomposition is not only the most 
natural way of connecting Ghh (|A4[) . to the radiative 
propagators g ra d ([!]) and g spp (O, as well as to recover 
previous results for the PEC Q, it is also related to the 
decomposition proposed in Ref. J25| in order to compare 
SPP with non-SPP mediated interaction. 

The real and imaginary parts of G^ d (R), G^(R), and 
G'spp(R) are compared with Ghh(R) hi Fig. [5] for the 
same parameters of Fig. |4j The most relevant feature 
observed in Figs. \5\ (a) and (b) is that the main con- 
tribution to Ghh(R) comes from the SPP propagator, 
G^ pp (R), which has a simple analytical form (|A9|) . This 



0,2 
0,1 

a -° 

a 

-0,1 



-0,2 



(a) 

A 


t 
1 
If 


Full 

— -RAD 

SPP 

„ EVAN 


■ 


f 1 

1 / 
I / 
/ / 

V / 

. i . 


\ 

i 


• 



0,5 1,0 1,5 2,0 2,5 3,0 3,5 4,0 

R/X 




0,5 1,0 1,5 2,0 2,5 3,0 3,5 4,0 

RIX 



FIG. 5. (Color). Real (a) and imaginary (b) parts of the 
propagator Ghh(R) (black solid line), as well as its con- 
stituent terms for radiative modes, G rad (-R) (red dashed line), 
SPP modes, Ggp P (R) (blue short-dashed line), and remaining 
evanescent modes G^(R) (green dash-dotted line). We use 
the same geometrical parameters of Fig. [4] 



allows us to find analytical expression for CI that will 
be presented in the next section. Notice that the agree- 
ment between Ghh(R) and G^ pp (R) has been previously 
reported for ID defects separated a distance larger that 
2-3A [13. Regarding non-SPP channels, G^ d (R) decays 
faster than G'l pp {R) being negligible small for R equal to 
a few X spp . On other hand, the real part of G^(R) is 
vanishing small (see Appendix), while its imaginary part 
is in anti-phase to G^ d (R). It must be noted that, as 
expected, the relative contribution of the different prop- 
agators changes when we approach the PEC limit [25H28J . 



A. Conditions for Interference 

In this section we compute the conditions for construc- 
tive and destructive interference of both out-of-planc 



and in-plane radiative powers. We start with the sim- 
pler of these two quantities, r\ ra d- Three simplifications 
help in finding the results for rj ra d- First, its interfer- 
ence pattern is accurately described by the normalized 
square field amplitude, rj ra d m |£^v| 2 , see Ufa). Second, 
Ghh{R) ~ G'sp p (R) 7 as we have learned from Fig. [3] 
Third, G^ pp (R) -C G s h- This last approximation is valid 
for R 3> X sp p, but we shall see it gives results that work 
surprisingly well even for R ~ X spp . Expanding E' N © 
into Mclaurin series oiG^ pp (R)/G s h and keeping only the 
leading term, we find 



E' 



N 



V 
-— xil-2Z E G spp (R), 



where Ze = E s h/I (|A7|) is the effective impedance of 
a single hole, which gives the modal amplitude at the 
hole opening as a function of the illumination. From 
the simplified expression for E' we can deduce that the 
interference pattern of the hole pair is set up by both the 
single hole impedance and the re-illumination of one hole 
by the other. 

The CI for i] ra d can be written in terms of the phase 
shift of both Z E and G™,(R). We thus define the 
single- hole phase shift, 4>ze, from Ze = \Ze\ ex p(*</ , z_E) J 
as well as the phase difference acquire by the SPP 
when traveling from one hole to the other, (f>hh, from 
G^pp(R) = \G^p P (R)\exp(i(j) hh ). An approximate ex- 
pression for the 4>hh can be obtained replacing the Han- 
kel function in G^ p (|A9[) by its asymptotic expression, 



H 



'(i) 



l 

hh 



(X) 



i (nx/2) 1 / 2 exp[i(x — tt/4)]. We have then 
,R — 7r/4. Keeping again the leading term in 



the expansion of \E' N \ , we obtain 



r/ rad sa 1 - 4 \Z E G^ p \ cos(k spp R + $ ZE - tt/4). (10) 

This equation clearly shows that the out-of-plane radi- 
ation depends both on the optical path traveled by the 
SPP when going to one hole to the other and the phase 
picked up by the field given the extra illumination pro- 
vided by the SPP coming from the other hole. The ap- 
proximate equation (|10j) is compared with full calcula- 
tions in Fig. UK a). We find that Eq. (JTDJl slightly under- 



estimates |£tv| 2 for R < X spp , but the agreement is excel- 
lent for R > X S pp- This nice agreement is related to the 
fact that non-SPP waves decays faster than SPP waves 
as a function of the distance, see Fig. [5J The leading role 
of SPP waves for large R have been already stressed in 
[251 - 4281 ] . Eq. (fl"0|) also agrees with the one proposed in Q 
following an intuitive interference plasmon model, which, 
in contrast to first-principles derivation of (|10[) , contains 
fitting parameters. 

It is straightforward to derive the CI of r\ ra d from Eq. 
(fTU|) assuming that the absolute value of G^ pp changes 
smoothly with i?, and that the dependence on R mainly 
comes from its phase. Then we have that extrema of r\ ra d 
appear at 



ftspp-H' 



•>ZE 



(11) 



where the integer value of n is equal to n = 2m — 1 for 
maxima, n = 2m for minima, and m = 1,2,3,...; while 
the condition for rj rac i = 1 is shifted in ir/2 with respect 
to the previous expression, i.e. 

7T 1 

k sppR- 7 + 0z_e = (m+ -)tt. (12) 

The single-hole phase shift, 0zb, is depicted in Fig. 
[6la) as function of the of the hole radius, r^, and for 
increasing metal thickness; h = 100 nm (blue dashed 
line), 150 nm (red solid line), and 250 nm (black short- 
dashed line). The large variation in 4>ze as function of 
Th (up to 7r/2 for increasing r% from 50 nm to 250 nm) 
accounts for the oscillations in r\ ra d observed in Fig. EHa). 
In Fig. Eta) we compare full calculations with the CI 
given by Eqs. (fTTj) and (fT2")l for the case r/, = 150, h = 
250 nm, for which <pzE = — 7r/4 (see Fig. |6ja)). An 
excellent agreement is obtained even for small values of 
R/X spp . Notice that r\ ra d is largely independent on the 
metal thickness h (although it is computed for a given 
value of of h in optically thick film) given that <f>zE is 
practically independent on h, see Fig. El 

Notice that the CI represented by Eqs. (fTTj) and (|T2l) . 
which are valid for a wide range of hole sizes (larger 
than the metal skin depth) and opaque metal films, are 
expressed in terms of the distance between the centers 
of the holes. A previous work [f| suggested that, for 
thin-metal films and small hole sizes, the CI are a func- 
tion of the edge-edge distance, independently from the 
hole radius. In our notation, this could only occurs if 
4>ze + 2rf l /\ spp — in Eq. (TT0|) . However, we observe in 
Fig. EJa) that —2rf l /\ spp (dash-dotted line) is equal to 
(/>ze only for a small region of the parameter space. This 
novel behavior demands further experimental work on 
opaque metal films and hole sizes larger than the metal 
skin depth. 

Similar CI can be developed for the in-plane scat- 
tered power rj S p P . As commented in the discussion of 
Fig. Bib), both terms l-E^I 2 and gf contribute to 
ri S pp in Eq. [5] We take the approximate expression of 
\E' N \ 2 from Eq. (jTUJ) and use the asymptotic expression 
gf = 1 + 2J[(k S ppR) found in the Appendix. We re- 
call that the last relation is exact when absorption is 
neglected, but otherwise it is still a good approximation. 
Using again the asymptotic expansion of the Bessel func- 
tion, the in-plane omittance is thus simplified to 



Vspp ~ I ^ 




Tt'^sppi*' 



cos(k spp R — 7r/4) 



.(13) 



This equation tells us that, given the normalized ampli- 
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FIG. 6. (Color), (a) Single-hole phase shift for the out-of- 
plane emittance, cj>ze, as function of the hole radius, rh, and 
for increasing metal thickness h; h — 100 nm (blue dashed 
line), 150 nm (red solid line), and 250 nm (black short-dashed 
line). The dash-dotted line represents the hole diameter nor- 
malized by X 3 pp- (b) SPP phase shift for the in-plane emit- 
tance of a single hole, <f>zp- 



order in G^L/Gsh, Tj sp p can be straightforwardly rewrit- 



ten to 



'Hspp 



1 



I y r~<hh 
\ Zj spp^' S pp 



cos(k spp R + 4> Z p - 7r/4),(14) 



where the effective impedance for the SPP channel 
Z spp — Ze — ^IGsppI) -1 takes into account both the 
excitation of the EM field inside the hole, characterized 
by Zei and the excitation of the SPP at the hole, given 
by (4|G*p„|) _1 . Like for the out-of-plane channel, the 



approximate Eq. (| 14[) shows an excellent agreement with 
full calculations in Fig. Sib). However, the behavior of 



tude of the electric field at the hole opening E' N , the (j)zp (defined from Z, 

interference pattern of the in-plane scattering power is 

determined by the SPP optical path between the two 

holes. However, in order to quantify the CI of r] S p P , we 

should include the modulation of the field given by I-E7J 2 

(flT)|) . Expanding the two terms in Eq. (jTB")) up to the first 



? spp = \Z spp \e^ zp ) differs from (j) Z E- 
Fig. EJb) renders §zp as a function of r/j, showing a char- 
acteristic peak centered near the cutoff radius r c = 168 
nm, cf. Fig. [5] The phase difference with respect to 
4>ze is about 7r/2 for r/j ^ r c , and decreases to zero for 
Th > r c . 



The extreme values of r] spp (p~4|) satisfy 
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r>'s 



}R ~4 



>ZP 



rnr, 



(15) 



where the integer value of n is equal to n — 2m for max- 
ima, n = 2m — 1 for minima, and m = 1,2,3, ... Notice 
that the values of rnr for rj spp are shifted in 7r with re- 
spect to the extreme values of r\ ra d (maxima are replaced 
by minima, and vise versa). This shifting is determined 
by the fact that the power traversing the hole is radiated 
into two complementary channels: r\ ra d and r\ spp . As for 
r]rad, the condition rj spp = 1 is shifted in 7r/2 with respect 
to the previous expression for extreme values, i.e. 



7T 1 

k sp pR - — + 4> Z p = (n + -)n. 



(16) 



In[3{b) we compare full calculations with the CI for the 
case rh = 150, h = 250 nm, for which (j>zp = 0.27 rad ps 
7r/4 (see Fig. [H[b)). As for cj>ze, an excellent agreement 
is obtained. 



V. CONCLUSIONS 

We have studied the emission pattern of a hole pair, 
focusing our attention in the role played by SPP reso- 
nances. Starting from the rigorous solution of the prob- 
lem, we have developed a SPP interference model that 
does not contain fitting parameters. This model provides 
simple analytical expressions for the interference pattern 
of both the out-of-plane and in-plane radiation channels, 
which nicely agree with full calculations for noble metals 
at optical frequencies. 

In agreement with experimental reports, both radiated 
powers oscillate with period \ spp - However, they show 
different trends as a function of the hole- hole distance and 
the hole radius. The amplitude of r\ ra d strongly oscillates 
with the hole radius, while the amplitude of r\ spp has a 
stronger dependence on R, but does not present such 
large variations with the hole size. 

Maxima of r) spp occur close to the conditions for con- 
structive interference of SPPs at the flat metal surface 
(R = m\ spp ) 1 while minima appear close to conditions 
for destructive interference of SPPs between the holes 
(R = (2m — l)X spp /2). CI for ri spp are shifted in it with 
respect to those of rj ra d (maxima are replaced by min- 
ima, and vise versa), because r\ ra d and r] spp are two com- 
plementary channels. The power traversing the hole is 
distributed into these two channels. 

We have also shown that two scattering mechanisms 
determine the interference pattern of the hole pair: (i) 
re-illumination by the in-plane SPP radiation and (ii) 
an effective impedance depending on the single-hole re- 
sponse. The conditions for interference only depend on 
the phase difference provided by each of the two scat- 
tering mechanisms. The large variation in the effective 
impedance of the single hole accounts for the oscillations 
of ri ra d as a function of the hole size. 
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Appendix A: Coupled-mode method 

In this section we briefly review the coupled-mode 
method for the optical transmission through holes, un- 
der the fundamental waveguide mode (TEn) approxi- 
mation. We refer to [lj| for the expressions of the full 
multimodc formalism and their derivation. Within the 
CMM, Maxwell's equations are solved self-consistently 
using a convenient representation for the EM fields. In 
both substrate and superstrate (see Fig. [l), the fields 
are expanded into an infinite set of plane waves with 
both p- and s-polarizations. Inside the holes the most 
natural basis is a set of circular waveguide modes [29|. 
The parallel components of the fields are matched at 
the metal/dielectric interface using surface impedance 
boundary conditions (SIBCs) [301. Although SIBCs ne- 
glect the tunneling of EM energy between the two metal 
surfaces, this effect is not relevant for a metal thickness 
larger than a few skin depths. 

At the lateral walls of the holes we choose the PEC 
approximation for the sake of analytical simplicity. We 
are thus neglecting absorption losses at the walls. Never- 
theless, we upgrade the PEC approximation introducing 
two phenomcnological corrections. First, the propaga- 
tion constant of the PEC fundamental mode is replaced 
by the one computed for a real metal. This improves the 
comparison between CMM and both experimental and 
FDTD results for both the spectral position of the peaks 
and the dependence of optical properties on the metal 
thickness. Second, enlarging the radius of the hole by 
one skin depth simulates the real penetration of in field 
into the metallic walls. This value for the enlargement 
provides the best agreement with FDTD simulations for 
an infinite periodic array of holes [21( . 

After matching the fields at the interface we arrive to 
the following system of tight binding-like equations 

[G sh - H]Ei{R) + G hh (R)E 2 (R) - G V E[(R) = h, 
[G sh - S] E[(R) + G hh (R)E' 2 (R) - G„E' 2 {R) = 0, 

where Ei is the modal amplitude of the electric field at 
the input opening of the i th hole, i = 1,2, and E[ is the 
same quantity but at the output opening. Two additional 
equivalent equations are needed for E 2 and E' 2 . Other 
relevant quantity is the illumination provided by the 
normal-incident p-polarized plane wave, with wavenum- 
ber k\ = 2-7r/A and admittance Y$ = y/e±, onto the 
lowest-energy mode 



h 



k> 



1 + z/ov^I' 



(Al) 



— 1/2 

where z s = e m is the metal impedance. In order to 
obtain a transmittance normalized by the flux impinging 
on the area covered by the holes, the illumination term 
/ already contains a factor (7rr^Y"o) -1 ^ 2 - The constant u 
satisfies J[ (u) = [29[ , where J(x) is the Bessel function 
of order 1, and the prime denotes derivation with respect 
to its argument 

The quantities E and G v represent scattering mecha- 
nisms already present in single holes. E is related to the 
bouncing back and forth of the waveguide fields inside 
the holes. Its value is 



E = Y„ 



#€!•*«* +/-e 



— Ci — ik z h 



li, 2 e< 



k z h 



Jw C 



ik z h 



(A2) 



where k z is the propagation constant of the waveguide 
mode, h is the metal thickness, /~ = 1 ± z s Y w , Y w = 
k z /k e2 is the admittance for the excited TEn mode, and 



e2k\. The quantity 



C_Tt/ ^-l 7/1 



/+ 2 e 



ik z h 



-&*& 



ik z h 



(A3) 



reflects the cou plin g between EM fields at the two sides 
of a given hole [19J . 

The propagator Ghh(R) represents the coupling of the 
two holes. It results from the projection of the Green's 
dyadic onto the waveguide modes in the holes. For the 
TEn mode, Ghh can be written as the following integral 
in the plane of the reciprocal space parallel to the metal 
surface 



Ghh(R) = Go 



G P (q) G s (q) 



qdq, (A4) 



where Go = 4fc 2 r 2 4 y / e/(u 2 — 1) and the two terms in the 
integrand represent the contribution of p- and s-polarized 
plane waves in the infinite semi-space in contact with 
the metal surface. The denominators of these two terms 
stand for the response of the metal plane. In particu- 
lar, the p-term has a pole at the SPP wavevector. The 
numerators G p and G s account for both the single hole 
response, which is a function of the hole radius ru, and 
the hole-hole interaction, a function of R. They read 



G p {q,r h ,R) 



Jf(k e qr h ) 
k^q 2 rf 



J[(hqR), 



G s {q,r hl R) 



k e qR 



(i-^y 



(A5) 
(A6) 



The integrand is written in adimcnsional units normal- 
izing the wavevector by k e = k\^/e. Notice that the R- 
dependent Bessel functions are obtained after the angu- 
lar integration in the kii = k e q(cos9, sin#) plane, where 
9 defines the direction of the component of wavevector 
parallel to the metal plane, kn . The dielectric constant e 
characterizes the dielectric material in contact with the 
metal surface (see Fig. [TJ e = e\ = 1 is used along this 
paper). 



The self-interaction term G s h is obtained after tacking 
the limit R — > in Ghh(R), i-e. using the identities 
lim x ^ J[{x) = linx^o Ji{x)/x = 1/2. 

Other relevant function is the effective impedance Ze, 
which is determined by the the three scattering mecha- 
nisms of the single hole (G s h, G„, and E), 






G s h 



(G s h 






(A7) 



fihh 
^rad 



We compute Ghh using the decomposition Ghh = 
G h ± + G™ ©, 



where G^ d 



spp 

contribution of radiative modes, G hh 



{R) represents the 
pp (R) designate the 



contribution of the plasmon pole to evanescent modes, 
and Gw(R) denotes the contribution of the remaining 
evanescent modes. The contribution of G^ d (R) can be 
written in terms of the functions gl1 d (R) and AG^ d (i?), 
which always take real values, 



G rad(R) = 9lld(R) + z ' s *AG Iad (R), 



where 



g^(R) = G / dq 



qq z Gp 

n -X- -/' 12 



qq z G s 



l 



q z z 



'12 



(A8) 



provides the interference term of the far field radiated 
from the holes ©, while the term proportional to the 
metal impedance z' s reads 



AG^ d (R) = Go f qdq 
Jo 



('» 



qlGs 



\q z + z' s \ 2 \l + q z z' s \ 



'12 



Both integrals are computed for free-propagating states 
(0 < q < 1). The real part of z' s is very small for typical 
noble metals, making Re[Gj^] 



I AG' 1 ' 1 



^ and Im[G^ dJ « 
z rad a good approximation for G 1£ [ d . The same re- 
lations hold for the single hole propagator, G s Jl d . 

For non-propagating states (q > 1) the integrand in 
Ghh(R) is prolonged into the complex q-plane, see [19| 
for details. The residue of the Cauchy integral gives the 
SPP wave confined to the metal/air interface 

G^. =mz' s Go J ' ikspprh) «' {1)l 



spp 



H'^'iksppR), (A9) 



spp h 



where k spp is the parallel component of the SPP wavevec- 
tor 



, ) » m (A) 
SPP ~ M e m (A)+e 



1/2 



(A10) 



Eq. (|A10|) defines k spp for a real metal and rigorous 
boundary condition at the metal/dielectric interface. In 
order to improve the accuracy our model we use Eq. 
(|A10|) instead of the approximate SPP wavevector for 
SIBCs, fc^ b p c = fc e [e(l - e^)] 1 / 2 . Besides the coupling 
propagator G^ pp , we define the radiative propagator g sp p, 
which provides the total SPP field radially scattered 
along all possible angular directions in the metal plane 



(J3j) . We have that g spp 
hole contribution read 



sh 
'spp 



g™' '• where the single- 



sh ._ 7rj4|aigo 

y spp c\ 



'Jl\fcspp r h) 



^spp^h 



and the interference term is equal to 



O*) 



Vspp 



(All) 



2Re [J[{k spp R)] + J{(2zlm[fc spp P]) - - 



2 
(A12) 



while 01 = |fc zp |Re[fc spp ]/(Im[fc zp ]|fc spp |), and k zp = (k 



k2 spp ) 1/2 - Notice that gf pp = g^/g%p © is independent 

Of TV 

For the sake of convenience, the integral for non-SPP 
evanescent states is computed along the vertical contour 
q = 1 + ih = <7_|_, /i£ [0, oo), after the integral variable is 
changed from q to h 



G 



hi. 



Go 
~2 



q + dh 



G ,, 



g: 



K z 
K Z G S 



K* - 

k*G! 



1 



IZ'K Z 



1 + iz'.K* 



(A13) 



where k, 



\/2i/i — /i 2 , and the Bessel function Ji(x) 
in both G p (|A5|) and G s (|A6|) is replaced by a Hankcl 

function of the first kind H[ (x). 

The propagator G%h is further simplified when the 
metal absorption is neglected, i.e. for Im[e m ] = Re[Zg] = 
and Im[Zg] = — \z'^. We find for the radiative modes 



that Rc[G^ d ] = g^ d , and Im[G&] 



\z 8 \AG% while 
for SPP modes Re[G spp (P)] = g S p P {R) = g%,[l + 
2J[(k spp R)]. For the remaining evanescent modes we ob- 
tain that Gpy is a pure imaginary function. The same 
relations that hold for Ghh are valid for G s h- Therefore, 
only the radiative and SPP terms contribute to the real 
part of in-plane propagator, 

Rc[G sh + G hh {R)\ = <? rad (P) + g spp (R). (A14) 

Under the lossless metal approximation, the total power 
traversing the two holes simplifies to Phoic = Re[G„PP'*] 
[19j | . We rewrite it in terms of E' with help of the relation 
G v E = {G sh + G hh (R)-T,}E',Le. 



P 



hole 



E'\ 2 Rc[G sh + G hh (R)]. 



As £ (IA2J1 is purely imaginary for a lossless media, this 
term does not contribute to Phoic- Using (|A14[) . we then 
have 



P 



hole 



P 



ad 



p.„ 



(A15) 



This equality represents the conservation of the power 
flux traversing the hole. These results can be easily 
generalized to an arbitrary number of holes, waveguide 
modes, and non-cylindrical geometries. It is also worth 
to mention that including absorption the computed pow- 
ers differ in less than 5% from the lossless case, even for 
a large number of defects [19l - |21] | . 

Finally, we recall that PEC is a particular case of a 
lossless metal with z' s = 0. In this case G^ d = O^d' w hile 
for non-propagating states of a PEC only G^ survives 



because G 



hh 
spp 
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